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ABSTRACT 

We present radial-velocity sequences acquired during three transits of the exoplanet HD 189733b and one transit of CoRoT-3b. We 
applied a combined Markov-chain Monte-Carlo analysis of spectroscopic and photometric data on these stars, to determine a full set 
of system parameters including the projected spin-orbit misalignment angle of HD 189733b to an unprecedented precision via the 
Rossiter-McLaughlin effect: f} = 0.S5°^^. This small but non-zero inclination of the planetary orbit is important to understand the 
origin of the system. On CoRoT-3b, results seem to point towards a non-zero inclination as well with ji = 37.6° ^?'3> t> ut tn i s remains 
marginal. Systematic effects due to non-gaussian cross-correlation functions appear to be the main cause of significant residuals that 
prevent an accurate determination of the projected stellar rotation velocity V sin(I) for both stars. 

Key words, binaries: eclipsing - planetary systems - stars: individual: HD 189733, CoRoT-3 - techniques: spectroscopic 



1. Introduction 

The spectroscopic transit, also known as the Rossiter- 



McLaughlin (RM) effect ( |Rossiter| p924| l ; |McLa ughlin ( 1924])), 
is a radial velocity (RV) anomaly superimposed on the radial- 
velocity curve arising from the Keplerian reflex orbit of the host 
star about its common centre of mass with a planet. As the planet 
transits, it covers - in the case of a prograde orbit - first the blue 
shifted part of the rotating star, shifting the overall spectrum 
slightly to the red. As the planet moves across the stellar disc, 
the radial velocity of the star's light centroid changes rapidly. 
The first RM effect caused by a transiting planet was observed on 
HD 209458 by |Queloz etal] pOOO) , In addition to the standard 
information that a photometric transit brings us, the RM effect 
permits us to measure the V sin{lj^of the star and is the only way 
to estimate the projected spin-orbit obliquity angle /3 (Gimenez 
2006b[ |Hosokawa|[T953| > (equal to -A ( |Ohta et al.||2005| ». This 
angle reflects the history of the planet and can therefore be used 
to constrain models of planetary orbital evolution. 



Send offprint requests to: Amaury. Triaud @unige.ch 

* using observations with the Harps spectrograph from the ESO 
3.6m installed at La Silla, Chile, under the allocated programmes 
072.C-0488(E) and 079.C-0828(A). The data is publicly available at 
the CDS (cdsarc .u-strasbg.fr) 

1 We used fhe |ohta et al.|(2005| notation differentiating the V sin(i) 
to the V sin(I). i is the projected inclination of the planet's orbit on the 
sky, whereas / is the projected inclination of the stellar equator on the 
sky. 



HD 189733 is a most interesting target due to its brightness 
and the possibility of observing it from both hemispheres fairly 
easily, while CoRoT-3b, orbiting a very fast rotating star and 
being itself in the middle of the Brown Dwarf desert (Deleuil 
et al.|2008| > (D08), is a unique object to study . Also HP 189733b 
has been extensively studie d (|Bouchy et al.|(|2005|l,|Winn et al. 
( 12006) 1 (W06), [Winn"eTaLl ( |2007| l (W07), |Boisse et al.| ( |2009> 
(B09)). The present paper is motivated by our recent acquisition 
of high-precision data - photometric for CoRoT-3, photometric 
and spectroscopic forHD 189733. Combining spectroscopic and 
photometric transits in one analysis allows us to refine the transit 
parameters, but also to ensure that the V sin(I) of the star and 
the spin-orbit angle /3, which are solely extracted from the RM 
effect, are tightly constrained. This is especially effective when 
analysing high precision data ; thus we use the high precision 
obtained on HD 189733 and present this star as a test-case for 
combining data sets. 

We describe the observations of these stars in section|2] then 
move to a description of the fitting process and its results with 
section[3]and explore a few reasons to explain the observed resid- 
uals that we obtained in section [4] Finally we will discuss the 
results and conclude. 



2. The Observations 

The observations were acquired using the High Resolution 
echelle spectrograph Harps mounted on the 3.6 m at the ESO 
observatory of La Silla, in Chile. The data were extracted using 
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the Data Reduction Software present at the telescope. A reanal- 
ysis of the data was performed later in Geneva, using the latest 
version of the software as in Mayor et al. (2009 1. 

One sequence of 1 1 Harps RV measurements in addition to 
those in the discovery paper (D08) was obtained on CoRoT-3 on 
August 26 th 2008 as part of the spectroscopic follow-up (072. C- 
0488(E)); all are around the transit. On average the new data has 
an estimated photon noise of 28.4m.* -1 . 

Four sequences were taken of HD 189733 on July 30 f/! , 
August 4" 1 and September % ,h 2006, and on August 29" 1 2007 
under the allocated programme 079.C-0828(A) and as part of the 
GTO, three of which are during transit. One sequence was taken 
off-transit to act as a comparison sequence. Out of the three RM 
sequences, two were obtained using a low cadence (one point ev- 
ery 10.5 minutes). One of these suffered from bad weather and 
one was taken with a high cadence (one point every 5.5 minutes). 
In total we have 78 new RV measurements including 37 during 
transit. The mean estimated uncertainty in the radial velocities 
due to photon noise is 0.98 m.s~ l . 



3. Fitting the Data 

Transiting planets have an important role to play in constraining 
planetary evolution models as well as atmosphere and interior 
models, therefore it is important that everything is done to en- 
sure that information extracted from the data is accurate. Both 
spectroscopic and photometric effects can be observed on a star 
experiencing a planetary transit. The two types of observation 
constrain parameters differently but arise from the same cause. 
Hence it is logical to fit both types of data simultaneously to de- 
termine a single set of parameters for the planet and ensure full 
consistency between the models. 



3.1. the modelling 

A code was developed using full Bayesian statistics in a Markov 
Chain Monte-Carlo (MCMC). The code is similar to the one de- 



scribe d in Collier |Cameron et al. ( 2007 1 and has already been 
used ( |Gillon et al.||2008[ |Bouchy et al.||2008] >. The philosophy 
here is to combine everything that is known about each star into 
the fitting process to better constrain the final result. So far, the 
periodic Doppler shift caused by the planet, the RM effect and 
a drift can be fitted to the spectroscopy. For the photometry, re- 
quiring data stripped of instrumental or stellar effects, primary 
and secondary transits can be fitted. Limb darkening coefficients 
can also be allowed to float if data of very high photometric pre- 
cision are available. 

The present version of the code fits up to 12 free parame- 
ters: the depth of the primary transit D, the RV semi-amplitude 
K, the impact parameter b, the transit width W, the period P, 
the middle of transit To, ecos(u>o) & e sin(ojo) (with e being the 
eccentricity and loq the angle between the line of sight and the 
periastron), V sin(I) cos(J3) & V sin(I) sin(J3), the RV drift F and 
the secondary transit depth Di- If some parameters are irrelevant, 
they can be fixed. These parameters have been chosen to reduce 
correlations and increase the exploration of parameter space. 

At each step i of the Markov chain, one set of j parameters is 
calculated from the previously accepted value (i — 1) following: 



P u = Pi-ij + f<r Pj G(0,l) 



(1) 



where Pj is a parameter, o>. is the 1 <x uncertainty, / is a factor 
ensuring 25% of steps are being accepted (see Tegmark et al. 
(2004)) and G(0, 1) is a random Gaussian number centred on 
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Fig. 1. top: Overall Keplerian fit of the RV data for CoRoT-3b. 
bottom: Composite plot showing both the spectroscopic and the 
photometric transit for CoRoT-3b. ( red) squares are Harps mea- 
surement which are not part of the RM sequence, (green) crosses 
are Sophie measurements, (blue) triangles show TLS data and 
(magenta) diamonds indicate the RM Harps sequence; (red) cir- 
cles are for the CoRoT photometry. 



zero with a standard deviation equal to 1. From these parameters 
a large variety of other - more physical - values can be inferred 
such as the stellar density p* (see Table[T|. 

The parameters are then used to calculate three different 
models. Both primary and secondary transit are calculated by 
using either codes developed by Man del & Agol| ( |2002| l, or by 
Gimenez (2006a i. The RV curve is modelled by the standard 
orbital equations ( |Mayor & Queloz|[l995| |Hilditch|[200T) and 
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Fig. 2. top: Overall Keplerian fit of the RV data for HD 189733b. 
bottom: RM effect on HD 189733. (red) squares shows Keck 
data from W06, (green) crosses are Sophie data from B09, (blue) 
triangles, (orange) circles, (magenta) diamonds are the new 
three Harps sequences on the RM effect and (black) oblique 
crosses show the off transit Harps sequence. Residuals are dis- 
played below and show a clear correlated signal during the tran- 
sit. 



by the code presented in Gimenez (2006b) for the RM effect. 
To model the stellar limb darkening, we use the quadratic law 
( Claret|2000 1. The model and the data are compared using a^ 2 
statistics. 

X 1 from the photometry is added to the value found for the 
spectroscopy. On that value Bayesian penalties are added. These 
can be estimated for every parameter for which we have indepen- 
dent prior knowledge of their value and error. The stellar mass 
M+ is allowed to vary subject to a Bayesian penalty on a value 
and its 1 cr error bar. All added, it creates a merit function: 



+ — -2 — + L 

M„ j=l 



(PiJ-Poj) 2 



(2) 



where here, Pj can be any parameter, fitted or physical and Pq 
is the value of the prior as M* is the floating parameter for stellar 
mass and M+ is the prior on the mass. 



This Qi is compared to the Qi-\, value calculated from the 
previous set of parameters, with the Metropolis-Hasting algo- 
rithm. This is repeated as many times as is necessary to ensure 
that the fitting has converged and that the exploration of param- 
eter space around the best value is truly random (meaning that 
the correlation length for each parameter is small compared to 
the number of accepted steps) and gives credible error bars. 

Obviously in addition to all these parameters, we also have 
to add one y velocity for each RV set and one normalisation 
constant for each photometric set; they are estimated by optimal 
averaging and optimal scaling in^ 2 calculations. The y velocity 
reflects the mean radial velocity due to the motion of the star in 
comparison to the Sun; its measured value varies with activity 
levels and between instruments by a small offset. The normal- 
isation factors re-estimate the normalisation of the lightcurves. 
The best fit parameters are chosen to be the set with the lowest 
X 2 and the error bars are calculated by taking the 68.3% lowest 
values of^- 2 and finding their extremes. 

3.2. the results 

One chain of 500,000 accepted steps was calculated for each star. 
Results for both stars appear in Table [T] The corresponding fits 
are displayed in Fig. [T] and Fig. [2] 

3.2.1. CoRoT-3b 

Nine free parameters, four y velocities and one photometric nor- 
malisation factor were used for CoRoT-3b to fit four RV se- 
quences totalling 29 measurements and one sequence of pho- 
tometry - the binned and phase folded data from D08 - with 400 
points in it. This amounts to 414 degrees of freedom. Results 
between this last paper and the present analysis are not very dif- 
ferent. The transit spectroscopic sequence covers little outside 
the RV anomaly of the RM effect and due to the faintness of the 
star and poor sampling during the transit, V sin(I) and /3 are not 
well defined. V sin(I) is found abnormally large at 35 ± 8 km.s~ l 
(see sect. |4.3[ >; /3 is different from zero only at the 2 cr level. 

Bayesian penalties were imposed for priors on the stellar 
mass M* = 1.37 ± 0.09 M G and also for the period P = 
4.25680 ± 0.000005 days because we made use of an already 
folded lightcurve as the data from D08 was released, rather than 
individual transits; the RM effect was not strong enough to con- 
strain P. To was allowed to float and is found to differ from the 
published value of 54283.1383+0.0003. This is possibly an arte- 
fact due to using a folded lightcurve. The RM sequence also suf- 
fers from a lack of continuum on either side of the transit. We 
employed the same limb darkening coefficients as D08. 

Our initial fit had a Bayesian penalty on the V sin(I) for the 
value of 17 ± 1 km.s~ l published in D08. Removing the penalty 
allowed us to find the current value of 35 + 8 km.s^ 1 andpermit- 
ted a minimisation of x 1 on the spectroscopy. See table^lwhere 
we show that having a V sin(I) = 35 km. s _1 and a [3 — 37.6° is 
a significant improvement of the model compared to an aligned 
system with the spectral analysis value of V sin(I) = 17 km. s~ l 
(because of a difference in P and To we also show results using 
parameters found by D08). 

We have 1 1 RV data points and three free parameters here 
(one y velocity, V sin(I) and /3) making the total number of de- 
grees of freedom eight. Thanks to the photometry, we have a se- 
cured detection of the planetary transit and of the Keplerian orbit 
and hence have a 100% chance that the RM effect will occur. It is 
only a matter of V sin(I) being different from zero. We obtain a 



4 



Triaud et al.: Rossiter-McLaughlin effect on CoRoT-3b & HD 189733b 



Table 1. Fitted and physical parameters found after fitting photo- 
metric and spectroscopic data of CoRoT-3b and of HD 189733b. 
The reduced^ 2 were calculated with 414 degrees of freedom for 
CoRoT-3b and 2890 for HD189773b. Nota Bene = -A; is 
used since the first reference of a projected spin-orbit angle was 
named thus, in|Hosokawa ( 1953 1. V sin(I)s are probably overes- 



timated, see section |4J for details 



Parameters (units) 



fitted parameters 



CoRoT-3b 



HD 189733b 



Table 2. Comparing x 1 f° r various solutions proposed for 
CoRoT-3b. 
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X re duced ~ ^ J ust f° r tne effect which has to be compared 
with &X 2 re duced ~ 5-962 if we adjust with a V sin(I) = 0. Hence, 
there is a clear detection of an RV deviation from the Keplerian 
orbit at the location of the RM effect. 



3.2.2. HD 189733b 

A total of 127 RV measurements were used for the fit includ- 
ing those published in B09 & W06 outside of the RM region. 
The RM effect in W06 is not used here. The RV data were fitted 
along with 2735 photometric points representing data from four 
transits in the z, (b+y)/2, V and / bands from W07. Compared to 
CoRoT-3b, 1 1 parameters were allowed to float, adding a drift F 
and P, plus two y velocities and three photometric normalisation 
factors more. This gives 2890 degrees of freedom. 

We applied a Bayesian penalty on M* alone, with the value 
0.82 + 0.03 M Q ( |Bouchy et al.|2005| l, i.e. using a & from eq. [2] 
with the two first terms only. We obtain very good constraints 
on P & Tq to a similar order of magnitude compared to Agol 



et al. ( 2008) > who used five Spitzer transits and four secondary 



transits. Most notable is the precision on the spin-orbit angle /3 = 
0.85° a 99.92% (3 cr) confident detection (W06 had found 
P — 1.4° ± 1.1). This value is found by comparing^ 2 with the 
value obtained by fixing /3 = and calculating the probability 
that indeed^ 2 was improved. 

HD 189733 is known to be active ( |Bouchy et aTJ ( |20"05] l), and 
stellar activity (stellar spots) causes changes in RV. Because the 
activity is acting on a timescale of the order of the rotational 
period, each of the four Harps sequences are expected to have 
no jitter related to stellar activity. But, the sequences were taken 
at different epochs, therefore at different activity levels. This is 
why each of the four sequences is fitted with its own y velocity. 
The change in y velocity of 33.5 m.s~ l is consistent with activity 
levels found by W06 and by B09. 

The results on HD 189733 confirm it is an exceptional tar- 
get against which to test the models, fitting techniques and data 
extraction. 



normalisation factors 
CoRoT 

FLWO 1.2mz-band 
T10 (b+y)/2 band 
MAGNUM 2 m V-band 
Wise 1 m 1-band 



0.9999998 



1.17 ±0.08 



0.99974 
0.99977 
0.99979 
0.99977 



1.21 ± 0.03 



4. Residual analysis 

It is clear when looking at the residuals of HD 189733b's RM ef- 
fect that a systematic error is present. The RMS of points within 
transit is 57 % larger than outside and their dispersion is corre- 
lated the same way for all three sequences. We also note that the 
residuals have a form very similar to those found by W06 (see 
Fig. 1 of that paper). 

On CoRoT-3 the RMS is comparable to the mean error bar 
outside and inside the transit indicating a good fit of the RV 
anomaly caused by the RM effect. 
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For both stars we find V sin(f) values larger than those 
present in the literature and found by an analysis of the spec- 
tral lines. For HD 189733, we have a 10% difference (or 2<x 
from the accepted value); for CoRoT-3, the value found is twice 
that which is inferred from spectral line analysis. We considered 
several possible causes of the effect observed in the residuals of 
HD 189733b and the large V sin(I) found on both stars. 

4. 1 . limb darkening 

We assumed that Harps was centred on the V band and limb 
darkening coefficients were chosen accordingly. The Harps 
spectral response was used to determine a new table of coeffi- 
cients in the manner described in Claret (2000). The difference 
between those newly found coefficients was small: u a = 0.6454 
instead of 0.6355 and u b = 0.1375 instead of 0.1488. for 
just the three HD 189773b RM Harps sequences, passed from 
213 to 206 ±20, consistent with no change at all. The shape in 
the residuals was not altered; this difference in limb darkening 
coefficients cannot explain the problem. 

Limb darkening being less well constrained near the limbs 
than at centre of the star, we omitted observations taken dur- 
ing ingress and egress, and re-fitted the RM effect. This tested 
whether the limb-darkening law itself could be causing the prob- 
lem. The fit was no better, the shape in the residuals was still 
present, ruling out this possibility as well. 

4.2. differential rotation 

The Sun exhibits surface differential rotation as a function of 
latitude. It is reasonable to assume that other solar like stars ro- 
tate differentially. The change in apparent rotational velocity for 
HD 189733b, thanks to a misaligned orbital angle and by tak- 
ing a value of differential rotation to be 10 % between the stellar 
poles and the equator, affects the amplitude of the RM effect very 
little: 



AVsin(I) = 0.3km.s~ 1 .AR il 



Sh |2. Jl-(b. cos P) 1 . tan/3 \ = 0.0218/?* 



6 V sin(I) = A V sin(I) .6h = 0.0065 km.s' 



(3) 



(4) 



(5) 



with h the altitude of the planet's path above the stellar equa- 
tor, projected onto the star (see Fig. [3]). The value is an order of 
magnitude lower than our error bars. 

It would also be expected to have an effect similar to a non 
zero spin-orbit angle: a difference in the amplitude on either side. 
Here the residuals are symmetric with respect to the centre of the 
transit. This said, with the level of precision now obtained on the 
spin-orbit angle of HD 189733b, it would be interesting to see 
how much of the /3 value is due to differential rotation, but this 
analysis goes beyond the scope of this paper. It can be empha- 
sised here that thanks to misaligned planets and in conjonction 
with more precise observations the Rossiter-McLaughlin effect 
will allow to detect differential rotation on stars without needing 
the spots that are used at the moment by Doppler tomography on 
stars like AB Dor. 

4.3. a systematic effect in the data reduction 




Fig. 3. Problem geometry and parameters used in equations 
present in section 



4.2 



Winn et aL) ( |2005| l (W05) pointed out that by fitting a symmet- 
ric Gaussian to a - by definition - non symmetric and varying 



cross correlation function (CCF) of the spectra in order to find 
the RVs and comparing these with a model which takes the cen- 
troid of the velocity weighted by the light emitted, we introduce 
a systematic error. The RV amplitude based on the symmetric 
Gaussian fit would be larger than that which a model with the 
same parameters would create. W05 tried to correct this effect 
by adding a polynomial on the RM model. This action was re- 
peated on W06 and subsequent papers. 

In W05, it was not shown what effect such a misuse of 
the model would have on the fit. Thus, we created two mod- 
els: one would act as the theoretical RM model does (Ohta's or 
Gimenez's), the other recreates the way the spectroscopic data is 
affected by the transiting planet, and the way the RVs composing 
the RM sequences are extracted. 

We created a grid on a star. Each element of the grid had two 
pieces of information on it: its intensity (taking limb darkening 
into account) and its apparent velocity on the line of sight. For 
the theoretical model, a planet was passed with similar charac- 
teristics as the two planets studied in this paper and the centroid 
of velocity weighted by the light was found for each position 
thereby recreating the RM effect. For the simulated data, a planet 
was also passed and a CCF of the star minus the contribution 
hidden by the planet was generated for each position assuming 
a Gaussian spectral line for each point on the star. A Gaussian 
function was fitted to that CCF and its minimum was taken as 
the simulated RV measurement. 

A comparison of the simulated data and the theoretical 
model showed what W05 had demonstrated: that the simulated 
data, using the same parameters, has a higher amplitude than 
the model predicts (see Fig. [4). Now, we want to see what a 
fit of that data would create. All transit parameters are heavily 
constrained by the photometry save the spin-orbit angle and the 
V sin(I). Assuming a /3 of zero we determined what would hap- 
pen if only V sin(I) was varying. 

On HD 189733, our test-case, by changing the V sin(I) value 
of the simulated data and comparing it to a model using the 
found value of 3.30 km. s^ 1 (see Table [TJ, we achieved a very 
good agreement when V sin(I) = 3.05km. s -1 . By subtracting 
one from the other, we obtained the theoretical residuals ex- 
pected for a fit such as that produced by the MCMC. The ob- 
served residuals from Fig. [2] were added on the graph (Fig. |4]i 
and showed a match, therefore indicating that a likely reason 
was found for these residuals. 
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Fig. 4. top: For HD 189733: dash (red) is the simulated spectro- 
scopic data with V sin(I) = 3.05 km. s~ l ; dash double dot (blue) 
is a RM model with V sin(I) = 3.30 km. dash dot (green) is 
the same model for a V sin(I) - 3.05 km. s , solid (magenta) is 
representing the residuals expected by subtracting the RM model 
by the simulated data; in (black) points, a sequence of residuals 
from Fig. [2] with Harps. These residuals are reproduced! bot- 
tom: For CoRoT-3b: dash (red) is the simulated spectroscopic 
data with V sin(I) = \1 km.s' 1 ; dash double dot (blue) is a RM 
model with V sin(I) = 35 km. s~ l ; dash dot (green) is the same 
model for a V sin(I) = 17 km.s , solid (magenta) is represent- 
ing the residuals expected by subtracting the RM model by the 
simulated data; with the (black) points, a sequence of residuals 
from Fig. [T]with Harps. 



The error bars presented in Table [T] seem to exclude such a 
large difference on V sin(I) at first sight, but it must be remem- 
bered that error bars depend on the model used; if the a pri- 
ori is wrong, then the a posteriori must be too. The presence of 
similar residuals in Fig. 1 of W06 shows that not everything is 
corrected with their method. Nevertheless, the value of V sin(I) 
(2.97 + 0.20 km. s^ 1 ) is similar to what this study infers showing 
the method developed by W06 is a good tool to estimate V sin(I). 



The data for HD 189733b were good enough to allow an 
adjustment. In the case of CoRoT-3b, where the sizes of error 
bars are large, such an adjustment was not possible. Instead we 
compared a theoretical model with the V sin(I) value found in 
Table[T]to simulated data with the value found by spectral analy- 
sis: 17+1 km.s~ l (D08). The expected residuals and the residuals 
of the Harps sequence were added on top (see Fig. Were it 
not for the abnormally large value of V sin(I) found by the fit, 
which gave clues that something might be going wrong, little 
information would have been extracted to support the idea the 
model is not adapted to the data because of error bars of the 
same order of magnitude as the RM effect itself, let alone the 
residuals! The absence of a clear signal due to a misfit of the 
model in the residuals is not enough to rule out the conclusions 
found with the case of HD 189733b that indeed the correlated 
residuals are caused by the imposition of a model not adapted 
to the data, or vice-versa that the data is not extracted with the 
same assumptions taken in the model. 

Two main paths are now being investigated to rectify the dis- 
covered problem: altering the models to take the effect into ac- 
count, which is the path already taken by W05 (but incomplete 
as the same residuals are observed), or finding alternative ways 
of looking at the data to extract it properly. It may be worth not- 
ing than tackling asymmetries in the CCF is underway in various 
areas, be it to understand stellar spots or classical cepheids stars 
(Nar detto et al.|2006) >. The RM effect is gaining popularity as a 
new powerful tool to estimate the rotation of stars. Potentially, 



on misaligned objects like XO-3b (Hebrard et al. 2008; Winn 



|et al.|2009| l HD 80606b (Mo utou et al.|2009| ) and CoRoT-3b, if 
confirmed, we could study differential rotation of stars other than 
the Sun. Before this can be achieved, however, we have to make 
sure that the RM effect is properly fitted. 



5. Discussion and Conclusions 

The results that we find are in accordance with previously pub- 
lished results, V sin(I)s aside. 

A spin-orbit angle is obtained for an exoplanet with unprece- 
dented precision with the value /3 - 0.85° and is a 3 <x detec- 
tion of an angle different from zero for HD 189733b's orbit. We 
confirm a marginal eccentricity which has been strengthened by 
the addition of the Harps data. It remains consistent with values 
found in |Agol et al. ( 2008b and B09. Both values are important 



to constrain planetary formation and evolution models. Results 
are consistent with the star having no radial velocity drift with 
time with F = 0.2+ 2 g m.s~ l .yr~ l . A comparison with a fit fix- 
ing T to zero showed that the variation in^- 2 was not significant: 
Ay 2 = 10 for an error on^ 2 of 30. The V sin(I) value found by 
the fit is spurious, as demonstrated in section |4~3| The real value 
is probably closer to 3.05 km.s^ 1 , value found for the simulated 
data to explain the theoretical models. 

On CoRoT-3b, the fit yields an implausibly high V sin(I) at 
35±8 km.s -1 , a two fold increase compared to the spectral analy- 
sis value, which is only explainable because of a discrepancy of 
assumptions between the model and the extraction of the data. 
T his overestimation of V sin(I) can also be read in the literature: 
in |Loeillet et~aT1 ( |2008| l a V sin(I) of 29.5 + 3 km.s- 1 is found 
for Hat-P-2 by fitting the RM effect without a prior. This value 
has to be compared to the value of 21 + 1 km.s- 1 found with an 
analysis of the spectral fines and of the photometry. Similarly for 
CoRoT-2b, in |Bouchy etaLl ([2008 ), a V sin(I) extracted from the 
RM model is found to be larger than its spectral and photomet- 
ric analyses counterparts. The discrepancy between the data and 
the model has for the moment mostly been observed in the case 
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of fast stellar rotators. Yet, this effect is also dependent on the 
Rp/R* ratio. This means that in the case where the RM effect is 
fitted independently of photometry leaving every parameter free, 
this mismatch between the centre of a fitted gaussian and the true 
velocity of the star's light centroid could lead to incorrect fitting 
of RM effects on slowly-rotating host stars of planets with high 
Rp/R+ ratios, and be mistaken for a V sin(I) problem. 

On CoRoT-3b we also get a marginal 97 % probability de- 
tection (equivalent to 2 cr) of an asymmetry in the current RM 
data. This result, if confirmed could shed a light into the debated 
origins and the uncertain nature of CoRoT-3b, a Brown Dwarf 
according to its mass and the presence of Deuterium burning, but 
a planet if created by accretion. A tormented history, reflected by 
a misaligned orbit, could point towards a more planetary origin 
for this body, notably via coalescence of large planetary bodies 
as outlined in Baraffe et al. (2008), who also point out the con- 
fusion of having Deuterium burning as the only limit between 
Brown Dwarfs and planets. The measure of /3 could become a 
way of segregating Brown Dwarfs from planets. 

Our study of HD 189733b began as a test-case for com- 
bining two different data sets. It also tested the RM model. 
CoRoT-3b came as an application of the analysis developed for 
HD 189733b and validates that a problem exists between the cur- 
rent models and the way the RV data is extracted. 

This paper shows the significance of fitting the various data 
sets in a combined way. In conclusion this combined analysis, 
(1) helps reduce the number of free parameters applied to fitting 
all the data, (2) breaks correlations between some parameters as 
they are fitted differently, hence exploring parameter space better 
and (3) insures a consistency on the transit parameters between 
the models and helps reveal model inconsistencies, as it did here 
with the case of the Rossiter-McLaughlin effect, therefore ensur- 
ing that systematics are not mistaken for physics. 
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